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Abstract 

We present the collective excitations spectrum (magnetoplasmon spectrum) of Dirac electrons in 
a weakly modulated single graphene layer in the presence of a uniform magnetic field. We consider 
electric modulation in one-dimension and the magnetic field applied perpendicular to graphene. We 
derive analytical results for the intra-Landau band plasmon spectrum within the self-consistent- 
field approach. We find Weiss oscillations in the magnetoplasmon spectrum which is the primary 
focus of this work. Results are presented for the intra-Landau band magnetoplasmon spectrum as 
a function of inverse magnetic field. These results are also compared with those of conventional 
2DEG. We have found that the Weiss oscillations in the magnetoplasmon spectrum are larger in 
amplitude compared to those in conventional 2DEG for the same modulation strength, period of 
modulation and electron density. 
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I. INTRODUCTION 



Recent progress in the experimental realization of a single layer of graphene has lead to 
extensive exploration of electronic properties in this system. Experimental and theoretical 
studies have shown that the nature of quasiparticles in these two-dimensional systems is 
very different from those of the conventional two-dimensional electron gas (2DEG) realized 
in semiconductor heterostructures. Graphene has a honeycomb lattice of carbon atoms. The 
quasiparticles in graphene have a band structure in which electron and hole bands touch at 
two points in the Brillouin zone. At these Dirac points the quasiparticles obey the massless 
Dirac equation. In other words, they behave as massless Dirac fermions leading to a linear 
dispersion relation ek = vk { with the characteristic velocity v ~ lO^m/s). This difference 
in the nature of the quasiparticles in graphene from conventional 2DEG has given rise to 
a host of new and unusual phenomena such as anamolous quantum Hall effects and a vr 
Berry phase These transport experiments have shown results in agreement with the 
presence of Dirac fermions. The 2D Dirac-like spectrum was confirmed recently by cyclotron 
resonance measurements and also by angle resolved photoelectron spectroscopy (ARPES) 
measurements in monolayer graphene Recent theoretical work on graphene multilayers 
has also shown the existance of Dirac electrons with a linear energy spectrum in monolayer 
graphene ^ . 

It was found years ago that if conventional 2DEG is subjected to artificially created 
periodic potentials in the submicrometer range it leads to the appearence of Weiss oscillations 
in the magnetoresistance. This type of electrical modulation of the 2D system can be carried 
out by depositing an array of parallel metallic strips on the surface or through two interfering 
laser beams . Weiss oscillations were found to be the result of commensurability of the 

electron cyclotron diameter at the Fermi energy and the period of the electric modulation. 

These oscillations were found to be periodic in the inverse magnetic field fl, y, 0. It 

is interesting to study the affect of electrical modulation of graphene layer on the Dirac 
electrons. Theoretical study of transport properties of Dirac electrons in graphene was 
recently carried out which showed the appearance of enhanced Weiss oscillations in the 
magnetoconductivity js]. 

One of the important electronic properties of a system are the collective excitations 
(plasmons). Weiss oscillations in the magnetoplasmon spectrum of conventional 2DEG has 
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been the subject of continued interest [9 
as in the eighties and more recently 



. Plasmons in graphene were studied as early 



111 ]. What distinguishes our work from earlier 



studies of plasmons in graphene is that we consider the effect of electric modulation on the 
magnetoplasmons which has not been studied before. In the present work, we study the 
affect of electrical modulation of a graphene monolayer on the magnetoplasmon spectrum 
within the self-consistent-field approach. We also analyze the dynamic, nonlocal dielectric 
function of this system and in particular consider the modulation-induced effects. 

The present paper is arranged as follows. In section II we give the formulation of the 
problem. We analytically determine the intra-Landau band plasmon spectrum in section 
III. In section IV we present the numerical results and discussion. Concluding remarks are 
made in section V. 



II. FORMULATION 

We consider two-dimensional Dirac electrons in graphene moving in the x-y-plane. The 
magnetic field {Bz) is applied along the z-direction perpendicular to the graphene plane. 
This system is subjected to weak electric modulation along the x-direction. We employ the 
Landau gauge and write the vector potential a.s A = (0, Bx, 0). The two-dimensional Dirac 



like Hamiltonian for single electron in the Landau gauge is = c = 1 in this work) 



Ho = va.{-iV + eA) (1) 

Here a = {cr^-, cr^jare the Pauli matrices and v characterizes the electron velocity. The 
complete Hamiltonian of our system is represented as 

H = Ho + U{x) (2) 

where Ho is the unmodulated Hamiltonian and U (x) is the one-dimensional periodic mod- 
ulation potential along the x-direction modelled as 

U{x) = VoCos{Kx) (3) 

here K = 271 /a, a is the period of modulation and Vq is the constant modulation amplitude. 
The Landau level energy eigenvalues without modulation are given by 

e{n) = v^/2eB \n\ = ujg^/\n\ (4) 



where ujg 



v\/2eBis the cyclotron frequency of the graphene electrons and n is an 



integer. The Landau level spectrum for Dirac electrons is significantly different from the 
spectrum for electrons in conventional 2DEG which is given as e{n) = uJc{n + 1/2), where 
ojc is the cyclotron frequency. 

Since we are considering weak modulation such that Vq is about an order of magnitude 
smaller than the Fermi energy Ep, we can apply standard perturbation theory to deter- 
mine the first order correction to the unmodulated energy eigenvalues in the presence of 
modulation with the result 



degeneracy of the Landau level spectrum of the unmodulated system with respect to ky 
is lifted in the presence of modulation. The formerly sharp Landau levels broaden into 
bands. The Landau bandwidths ~ oscillates as a function of n since Ln{u) are oscillatory 
functions of the index n. Similar features were noted in the Landau level spectrum of 
conventional electrons in density modulated 2DEG. The difference in the case of density 
modulated Dirac electrons is firstly that the Landau level spectrum depends on the square 
root of both the magnetic field B and the Landau band index n as against linear dependence 
in case of conventional electrons. Secondly, for Dirac electrons the energy eigenvalues in the 
presence of modulation given by Eq. (6) contain a term which is a linear combination of 
Laguerre polynomials with indices n and n — 1 while for conventional electrons there is only a 
single term that contains Laguerre polynomial with index n. We expect that this modulation 
induced change in the electronic density of states to influence the dielectric response and 
the collective excitations of the graphene system. 

III. INTRA-LANDAU-BAND PLASMON SPECTRUM OF GRAPHENE 




(5) 





The dynamic and static response properties of an electron system are all embodied in 
the structure of the density- density correlation function. We employ the Ehrenreich- Cohen 
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self-consistent-field (SCF) approach |12[| to calculate the density-density correlation function. 
The SCF treatment presented here is by its nature a high density approximation which has 
been successfully employed in the study of collective excitations in low-dimensional systems 
both with and without an applied magnetic field. It has been found that SCF predictions 
of plasmon spectra are in excellent agreement with experimental results. Following the SCF 
approach, one can obtain the dielectric function as 

e(g, u) = l- v^{q)U{q, u). (7) 

where Vc{q) = Q = {qI + QyY^'^, is the background dielectric constant and Il{q,uj) is 



Kq 

the non-interacting density-density correlation function 







X [e{n,Xo + x'q) - 6{n',Xo) + uj + ir]] , (8) 
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where Xq = l^ky, x'q = l^qy, C„„/(^) = (^2!/^!!) i-^) "^1 with rii 

max{n,n'),n2 = min(n, n') and L'^{x) an associated Laguerre polynomial. In writing the 
above equation we converted the ky-sum into an integral over xq by the following relation: 

a 

-jY^ 2^;^ J dxQ ■ ■ ■ . This is a convenient form of n(g, c<j) that facilitates writing of 

ky 
the real and imaginary parts of the correlation function. 

The real and imaginary parts of the density- density correlation function are the essential 

ingredients in the theoretical considerations of such diverse problems as high frequency 

and steady state transport, static and dynamic screening and correlation phenomena. The 

plasmon modes are determined from the roots of the longitudinal dispersion relation 

l-t;,(g)Ren(g,cu) = 0, (9) 

along with the condition Imn(g, ct;) = to ensure long-lived excitations. Employing Eq.(8), 
Eq.(9) can be expressed as 

n,n' 

with 

/i(-) = p hx„- — „ (11) 

J e{n,xo) - e{n,xo + x'q) + uj 



{P denotes principal value integral). The plasmon modes originate from two kinds of elec- 
tronic transitions: those involving different Landau bands (inter-Landau band plasmons) 
and those within a single Landau-band (intra-Landau band plasmons). 

Inter-Landau band plasmons involve the local ID magnetoplasma mode and the 
Bernstein- like plasma resonances [9], all of which involve excitation frequencies greater than 
the Landau-band separation. On the other hand, intra-Landau band plasmons resonate 
at frequencies comparable to the bandwidths, and the existence of these modes is due to 
finite width of the Landau levels. The occurrence of such intra-Landau band plasmons is 
accompanied by SdH type of oscillatory behavior in 1/B. These SdH oscillations occur as a 
result of emptying out of electrons from successive Landau bands when they pass through 
the Fermi level as the magnetic field is increased. The amplitude of these oscillations is a 
monotonic function of the magnetic field, when the Landau bandwidth is independent of the 
band index n. In the density modulated case, the Landau bandwidths oscillate as a function 
of the band index n, so it is fair to expect that such oscillating bandwidths would affect 
the plasmon spectrum of the intra-Landau band type, resulting in new kind of oscillations 
called Weiss oscillations studied extensively in conventional 2DEG. Here, we study these 
Weiss oscillations in the magnctoplasmon spectrum in graphene. 

To examine the intra-Landau-band plasmon spectrum we consider electronic transitions 
within a single Landau band. In this case <^ £„, this approximation enables us to expand 
the argument of the distribution function about and substituting this expansion into 
Eq.(ll) yields 

a 

/i(cu) = " y dxof{en + Fn cos{Kxq)) 



X (1 - —{cos{Kxo) - cos{K{xo + x'o))]). (12) 

The first integral in the above equation vanishes while we fold the range of integration from 
to a/2 by multiplying a factor of 2 in the second integral. Hence wc can express the 
intra-Landau band plasmon dispersion relation, from Eqs. (10, 11, 12) and C„„/(a;) — > 1 for 
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n — n , as 



16e2 



a/2 



sm (-Xo) 



X 



y^Fn J dxof{en + Fn cos{Kxq)) cos{Kxq) 



(13) 



This can be written as 



2 16e2 1 . 2/7r , 

a; = — -^sm {-Xq)x 

kq nar a 



Y.\Fn\A„ 



(14) 



here Ar 



a \Pn 
2-K F„ 



a/1 — ^^(1 — A„) which results from the analytic evaluation of the Fermi 
function integral in Eq.(13) that can be carried out at T = O.The integral contributes only 



with Ep the Fermi energy. Hence the zero 



for A„ < 1, where A„ is given by A, 
temperature intra-Landau-band plasmon dispersion relation reduces to 

o 8e2 1 



kq nal'- 



sm 



;-x(,)x5]|F„iyr^^(i-A„). 



(15) 



In general, the inter- and intra-Landau-band modes are coupled for arbitrary magnetic field 
strength. So long as < vy/neB, mixing of the inter-and intra- Landau band modes 
is small. Only the intra-Landau-band mode (cu) will be excited in the frequency regime 
vy/neB > uj |F„| . 



IV. NUMERICAL RESULTS AND DISCUSSION 

The intra-Landau-band plasma frequency for Dirac electrons given by equation (15) 
is shown graphically (solid curve) in Fig.l as a function of 1/B. In the same figure we 
also display the intra-Landau-band plasma frequency for standard electrons in conventional 
2DEG (dashed curve). The following parameters were employed for graphene: k — 2.5, 

= 3 X 10^^ m-^ v = 2.6 eV A, a = 382 nm and Vq = 0.5 meV. We also take g^; = and 
qy = .Olkp, with kp = {txudY^'^- For conventional 2DEG (a 2DEG at the GaAs-AlGaAs het- 
erojunction) we use the following parameters: m = .07me(me is the electron mass), k = 12 
and Ud = 3 x 10^^ m~^ with the modulation strength and period same as in the graphene 
system. In Fig.l, the modulation induced oscillations due to the ID modulation are clearly 
visible superposed on the sharp SdH-type oscillations. The intra-Landau-band plasmons 
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have frequencies comparable to the bandwidth and they occur as a result of broadening of 
the Landau levels due to the modulation in our system. These type of intra-Landau-band 
plasmons accompanied by regular oscillatory behavior (in 1/B) of the SdH type was first 
predicted in [13] for tunneling planar superlattice where the overlap of electron wavefunction 
in adjacent quantum wells provides the mechanism for broadening of Landau levels. The 
SdH oscillations occur as a result of emptying out of electrons from successive Landau levels 
when they pass through the Fermi level as the magnetic field is increased. The amplitude of 
these oscillations is a monotonic function of the magnetic field when the Landau bandwidth 
is independent of the band index n. In the density modulated case, the Landau bandwidths 
oscillate as a function of the band index n, with the result that in the plasmon spectrum 
of the intra-Landau band type, there is a new kind of oscillation called Weiss oscillation 
which is also periodic in but with a different period and amplitude from the SdH type 
oscillation. 

The origin of these two types of oscillations can also be understood by a closer analytic 



examination of equation (15). In the regime, vyneB > the unit step function vanishes 
for all but the highest occupied Landau band, corresponding, say, to the band index A^. 
The sum over n is trivial, and plasma frequency is given as ~ |Fjv|"'^^^ (1 — A'j^y^^6(l — 
Ajv). The analytic structure primarily responsible for the SdH oscillations is the function 
6(1 — An), which jumps periodically from zero (when the Fermi level is above the highest 
occupied Landau band) to unity (when the Fermi level is contained with in the highest 
occupied Landau band). On the other hand, Weiss oscillations which represent the periodic 
modulation of the amplitude of SdH oscillation are largely a consequence of the oscillatory 
factor |F/v|^''^. 

From the results displayed in Fig.l, we can also make a comparison of the Weiss os- 
cillations in the magnetoplasmon spectrum of Dirac electrons in graphene and standard 
electrons in conventional 2DEG. For the same modulation strength, period of modulation 
and electron density we find that the Weiss oscillations in graphene have a larger amplitude 
compared to conventional 2DEG. This can be attributed to the larger characteristic velocity 
(f ~ lO^m/s) of electrons in graphene compared to the Fermi velocity of standard electrons 
and smaller background dielectric constant k in graphene in contrast to conventional 2DEG. 

The SCF theory presented here is expected to be a good quantitative approximation 
since the dimensionless electron density parameter in graphene is a constant (< 1) 11| . 
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hence we are in the high density regime. It should also be noted that for electron densities 
under consideration the effective Fermi temperature Tp = Ep/kB in graphene is very high 
(~ ISOOX for ?T,£) ~ 10^^ ni"^), so our T = theory is expected to be valid at higher 
temperatures. The only constraint is that the intra-Landau band plasmons are due to 
electronic transitions within a single Landau band and the existence of these modes is due 
to finite width of the Landau levels caused by the modulation. Therefore modulation induced 
effects will be observable for temperatures where the thermal energy is less than the Landau 
band width such that ksT < 2|F„|. In effect the theory is valid for experimental situations 
where ksT is smaller than ujg, Fn and ep- 

V. CONCLUSIONS 

We have determined the intra-Landau band plasma frequency for Dirac electrons in a 
modulated graphene layer in the presence of a magnetic field employing the SCF approach. 
We show that the intra-Landau band magnetoplasmon spectrum exhibits the electronic 
density of states modulation in the form of oscillating magnetoplasmon frequencies (Weiss 
oscillations). We also compare these results with those of the conventional 2DEG. We have 
found that the Weiss oscillations in the magnetoplasmon spectrum arc larger in amplitude 
compared to those in conventional 2DEG for the same modulation strength, period of mod- 
ulation and electron density. 
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